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We investigate the chiral condensate and the dressed Polyakov loop or dual chiral condensate 
at finite temperature and density in two-flavor Nambu-Jona-Lasinio model. The dressed Polyakov 
loop is regarded as an equivalent order parameter of deconfinement phase transition in a confining 
theory. We find the behavior of dressed Polyakov loop in absence of any confinement mechanism 
is quite interesting, with only quark degrees of freedom present, it still shows an order parameter 
like behavior. It is found that in the chiral limit, the critical temperature for chiral phase transition 
coincides with that of the dressed Polyakov loop in the whole (T, fi) plane. In the case of explicit 
chiral symmetry breaking, it is found that the transition temperature for chiral restoration T c x is 
smaller than that of the dressed Polyakov loop Tf in the low baryon density region where the 
transition is a crossover. With the increase of current quark mass the difference between the two 
transition temperatures is found to be increasing. However, the two critical temperatures coincide 
in the high baryon density region where the phase transition is of first order. We give an explanation 
on the feature of T* = T,? in the case of 1st and 2nd order phase transitions, and T* < T^ in the 
case of crossover, and expect this feature is general and can be extended to full QCD theory. Our 
result might indicate that in the case of crossover, there exists a small region where chiral symmetry 
is restored but the color degrees of freedom are still confined. 
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I. INTRODUCTION 

Quantum Chromodynamics (QCD) in nonperturbative 
regime has two fundamental properties: chiral symme- 
try is dynamically broken and color charges are con- 
fined. The interplay between chiral symmetry breaking 
and confinement as well as the chiral and deconfinement 
phase transitions at finite temperature and density are 
of continuous interests 1-6]. These two transitions are 
characterized by the breaking and restoration of chiral 
and center symmetry, which are well defined in two ex- 
treme quark mass limits. In the chiral limit when the 
current quark mass is zero m — 0, the chiral condensate 
is the order parameter for the chiral phase transition. 
When the current quark mass goes to infinity m — > oo, 
QCD becomes pure gauge SU(3) theory, which is cen- 
ter symmetric in the vacuum, and the usually used order 
parameter is the Polyakov loop [l| . 

The relation between the chiral and deconfinement 
phase transitions is very important for the phase diagram 
at high baryon density region [7] . It is conjectured in Ref. 
8] that in large N c limit, a confined but chiral symmetric 
phase, which is called quarkyonic phase can exist in the 
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high baryon density region. It is very interesting to study 
whether this quarkyonic phase can survive in real QCD 
phase diagram, and how it competes with nuclear matter 
and the color superconducting phase. (It is noticed that 
in Ref. 9], it is found that at zero chemical potential, the 
lattice results for the thermodynamical properties have a 
very mild dependence on the number of colors.) 

Lattice QCD at the current stage cannot go to very 
high baryon density. For zero chemical potential, pre- 
vious lattice results show that these two transitions oc- 
cur at the same temperature, e.g, in Ref. [HM141]. and 
also in review papers [la, [pq|. In recent years, three 
lattice groups, MILC group [17j . RBC-Bielefeld group 
[lc| and Wuppertal-Budapest group [19rf21| have stud- 
ied the chiral and deconfinement phase transition tem- 
peratures with almost physical quark masses. The RBC- 
Bielefeld group claimed that the two critical tempera- 
tures for N f = 2 + 1 coincide at T c = 192(7) (4)MeV. 
The Wuppetal-Budapest group found that for the case 
of Nf = 2 + 1, there are three transition tempera- 
tures, the transition temperature for chiral restoration 
of u,d quarks T c x( " d) = 151(3) (3) MeV, the transition 
temperature for chiral restoration of s quark T c x = 
175(2) (4)MeV and the deconfinement transition temper- 
ature T c d = 176(3) (4)MeV. From this result, we can read 
that the critical temperatures for different transitions are 
different. According to the Wuppetal-Budapest group, 
this is the consequence of the crossover nature. 

In the framework of QCD effective models, there is still 



no dynamical model which can describe the chiral sym- 
metry breaking and confinement simultaneously. The 
main difficulty of effective QCD model to include con- 
finement mechanism lies in that it is difficult to calcu- 
late the Polyakov loop analytically. Currently, the pop- 
ular models used to investigate the chiral and deconfinc- 
ment phase transitions are the Polyakov Nambu-Jona- 
Lasinio model (PNJ L - ) [22H 27J and Polyakov linear sigma 
model (PLSM) [2|| [29|. However, the shortcoming of 
these models is that the temperature dependence of the 
Polyakov-loop potential is put in by hand from lattice 
result, which cannot be self-consistently extended to fi- 
nite baryon density. Recently, efforts have been made 
in Ref.|30] to derive a low-energy effective theory for 
confinement-deconfincment and chiral-symmetry break- 
ing/restoration from first-principle. 

Recent investigation revealed that quark propagator, 
heat kernels can also act as an order parameter as they 
transform non trivially under the center transformation 
related to deconfinement transition |3ll433| . But the ex- 
citing result is the behavior of spectral sum of the Dirac 
operator under center transformation [32], |3J-|36[ . A new 
order parameter, called dressed Polyakov loop has been 
defined which can be represented as a spectral sum of the 
Dirac operator [36|j . It has been found the infrared part of 
the spectrum particularly plays a leading role in confine- 
ment [32J . This result is encouraging since it gives a hope 
to relate the chiral phase transition with the confinement- 
deconfinement phase transition. The order parameter for 
chiral phase transition is related to the spectral density 
of the Dirac operator through Banks-Casher relation [4j. 
Therefore, both the dressed Polyakov loop and the chiral 
condensate are related to the spectral sum of the Dirac 
operator. 

Behavior of the dressed Polyakov loop is mainly stud- 
ied in the framework of Lattice gauge theory [37|, |38| . 
Apart from that, studies based on Dyson- Schwinger 
equations |39J - l4l| and PNJL model [42| have been car- 
ried out. In those studies the role of dressed Polyakov 
loop as an order parameter is discussed at zero chemical 
potential. In this work, we will investigate the decon- 
finement phase transition by using the dressed Polyakov 
loop as an equivalent order parameter in the Nambu- 
Jona-Lasinio (NJL) model. As we know, the NJL model 
lacks of confinement and the gluon dynamics is encoded 
in a static coupling constant for four point contact inter- 
action. However, assuming that we can read the infor- 
mation of confinement from the dual chiral condensate, 
it would be interesting to see the behavior of the dressed 
Polyakov loop in a scenario without any explicit mecha- 
nism for confinement. 

In this work we study the phase transitions in the NJL 
model in (T, /i) plane in chiral limit as well for small 
quark mass limit. This paper is organized as follows: We 
introduce the dressed Polyakov loop as an equivalent or- 
der parameter of confinement-deconfinement phase tran- 
sition and the NJL model in Sec. |TTJ Then in Sec JIIH 
we investigate the phase transitions in T — fj, plane in 



the chiral limit and in the case of explicit chiral sym- 
metry breaking, respectively. We offer an analysis on 
the relation between the chiral and deconfinement phase 
transitions in Sec. IIV1 At the end, we give the conclusion 
and discussion. 



II. DRESSED POLYAKOV LOOP AND THE NJL 
MODEL 

Before going to our work, let us briefly introduce the 
dressed Polyakov loop. To do this we have to consider a 
U(l) valued boundary condition for the fcrmionic fields 
in the temporal direction instead of the canonical choice 
of anti-periodic boundary condition, 



i>(x,p) = e- i *i>(x,0), 



(1) 



where < <j) < 2-7T is the phase angle and ft is the inverse 
temperature. 

Dual quark condensate S n is then defined by the 
Fourier transform (w.r.t the phase <f) of the general 
boundary condition dependent quark condensate [36l43a | , 



Sn = " / 7T- 
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W)< 



(2) 



where n is the winding number. 

Particular case of n = 1 is called the dressed Polyakov 
loop which transforms in the same way as the conven- 
tional thin Polyakov loop under the center symmetry and 
hence is an order parameter for the deconfinement tran- 
sition |36l - [38| . It reduces to the thin Polyakov loop and 
to the dual of the conventional chiral condensate in infi- 
nite and zero quark mass limits respectively, i.e., in the 
chiral limit m — - > we get the dual of the conventional 
chiral condensate and in the m — >■ oo limit we have thin 
Polyakov loop |36l - l38j . 

The two-flavor NJL model is described by the La- 
grangian density in the form of 



C = ^{i'fd^. - m)^ + G s [(^) 2 + (^75T^) 2 



(3) 



where only scalar and pseudo-scalar channels are con- 
sidered, to is the current quark mass. The (f> dependent 
thermodynamic potential in the mean field level for scalar 
four point interaction is given by, 



n 



NJL 



-2N 



fl A ^[^ P + lln(l + e-^-) (4) 



(2vr) 3 L "" p '' ft' 



3 



hIn(l + e-^)]+G s ( ff )i 



Where, Nf is the number of flavors which is 2, E p = 
E p ±fi±T&, (a)^ = (V^VV- The single particle energy is 



&P = V I V I + M 2 an d constituent mass is given by M = 
777 — 2G s (a) ( j } . The last term T$ in the expression for 
E p arises because of the imposed U{\) valued temporal 



boundary condition parameterized by the angle (f> and is 
given by $ = — in + i(f>. 

We work in the isospin symmetric limit and contribu- 
tion from pions arising out of pseudo scalar interaction 
is not taken into account. The thermodynamic potential 
contains imaginary part. We take only the real part of 
the potential and the imaginary phase factor is not con- 
sidered in this work. The mean field (cr)^ is obtained 
by minimizing the potential for each value of </> S [0, 2ir) 
for fixed T and /i. The conventional chiral condensate is 
(a)„ = (V^tt- The dressed Polyakov loop Si is obtained 
by integrating over the angle. The values of the parame- 
ters A and G s are taken as 0.6315GeV and 5.498GeV~ 2 , 
respectively. 
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III. NUMERICAL RESULTS 
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FIG. 1: Angle variation of (a)^, for different values of temper- 
atures and chemical potentials in the case of chiral limit. The 
solid line corresponds to T — 150MeV,/i = lOOMeV, dashed 
line corresponds to T = 250MeV,/i = lOOMeV, dash-dotted 
line corresponds to T = 40MeV, [i = 300MeV, and dot corre- 
sponds to T = 150MeV,At = 300MeV. 

In this work, we investigate phase transitions for two 
cases, i.e., in the chiral limit and with explicit chiral sym- 
metry breaking. Fig. Q] and [2] show the behavior of the 
angle dependence of the general chiral condensate for var- 
ious chemical potentials and temperatures for to = and 
to = 5.5MeV, respectively. The four curves presented 
in each figure represent two temperatures above and be- 
low the critical temperature for two particular values of 
the chemical potential. Same qualitative features have 
been found for both the quark masses. The variation is 
symmetrical around <f> — -k as reported in other studies 

[MSI- 

Almost no variation with respect to angle is found fol- 
low temperatures. As the temperature increases the vari- 
ation over the angle grows. We expect the absolute value 
of the chiral condensate decreases with the increase of 
temperature. However, from the figure, this conventional 



FIG. 2: Angle variation of {<?},/, for different values of tempera- 
tures and chemical potentials in the case of m — 5.5MeV. The 
solid line corresponds to T = 150MeV,/i = lOOMeV, dashed 
line corresponds to T = 250MeV,/i = lOOMeV, dash-dotted 
line corresponds to T = 20MeV, fj, = 340MeV and dotted line 
corresponds to T - 40MeV, fi - 340MeV. 



behavior of the chiral condensate with temperature only 
persists up to a certain angle, beyond which the oppo- 
site behavior is observed. The plateau around <f> = tt is 
more flat above T c in case of zero current quark mass and 
its value is consistent with the expectation of complete 
restoration of chiral symmetry in the chiral limit. 
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FIG. 3: The conventional chiral condensate ~(<j)tt and the 
dressed Polyakov loop Si as functions of temperature for dif- 
ferent values of the chemical potentials in the chiral limit. 
From right to left the values of the chemical potential are 
0, 200, 300MeV, respectively. 

Fig. [3] and S] show the behavior of the conventional chi- 
ral condensate —(<t)tt and the dressed Polyakov loop Ei 
at different chemical potentials as functions of tempera- 
ture for to = and m = 5.5MeV, respectively. For both 
cases, it is observed there are three temperature regions 
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FIG. 4: The conventional chiral condensate — (cr) w and the 
dressed Polyakov loop Ei as functions of temperature for dif- 
ferent values of the chemical potentials in the case of explicit 
chiral symmetry breaking m = 5.5MeV. From right to left 
the values of the chemical potential are 0, 200, 300MeV, re- 
spectively. 



for —{<j)n and Si. For —(a) n ^ at smaller temperatures it 
remains constant at a value corresponding to the value of 
the conventional chiral condensate in the vacuum, then 
it rapidly decreases in a small window of temperature 
and eventually almost saturates to a lower value. The 
rapid decreasing occurs at different temperatures for dif- 
ferent values of the chemical potentials. On the other 
hand the behavior for the dressed Polyakov loop is just 
the opposite. It remains almost zero for small temper- 
atures and then rises rapidly, finally saturates to a high 
value which varies very slowly with temperatures. The 
almost zero value of Si for small temperatures is due 
to the fact that the U(l) boundary condition dependent 
general quark condensate nearly does not vary with the 
angle <fi for small temperatures (see Eq. [2J . 

For finite quark mass, near the critical temperature 
region, both — (cr)^ and Si change more slowly than those 
in the case of chiral limit. 

Fig. [5]and[6]show the T — fi phase diagram for the case 
of to = and to — 5.5MeV, respectively. The transition 
temperatures are calculated from the slope analysis of 
the conventional chiral condensate (er)^ and the dressed 
Polyakov loop. The transition temperatures calculated 
from the conventional chiral condensate represent the 
chiral phase transition temperature. On the other hand 
the behavior of the dressed Polyakov loop is supposed 
to indicate the deconfinement transition temperature. In 
our present framework confinement is not accounted for, 
however, if we look at the curves presented in figure [3] 
and|U they still show an order parameter like behavior. 

For to = case, we find almost exact matching for the 
transition temperatures calculated from these two quan- 
tities in the whole T — /j, plane as shown in Fig. [5] 

For the case of finite quark mass m = 5.5MeV, it is ob- 
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FIG. 5: T — fi phase diagram for the case of chiral limit. The 
solid line is the critical line for Ei, and the dashed line is the 
critical line for conventional chiral phase transition. The solid 
circle indicates the critical point for chiral phase transition. 
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FIG. 6: T — n phase diagram for the case of m = 5.5MeV. 
The solid line is the critical line for Ei , and the dashed line is 
the critical line for conventional chiral phase transition. The 
solid circle indicates the critical end point for chiral phase 
transition. 



served from Fig. [5] that the two critical temperatures are 
different in the low baryon density region. The difference 
however decreases from low to high chemical potential, 
and the two critical temperatures start to match around 
the critical end point for chiral phase transition. For 
zero chemical potential and small current quark mass 
to = 5.5MeV, we find about 7MeV difference between 
T* and Tf, and Tf < Tf '. Similar trend has been ob- 
served in another study based on Dyson-Schwinger ap- 
proach |39|, where they found chiral transition to occur 
about 10 — 20MeV below the deconfinement transition. 
Though these studies are not a complete one and these 
difference may be due to the effects of crossover transi- 



tion. As pointed out in [21| during crossover, different 
observables are expected to behave differently and there 
is no way to define a unique crossover temperature. 

We extend our study further to see what happens if 
we increase the current quark mass further. We find at 
zero chemical potential, the difference between the transi- 
tion temperatures calculated from dressed Polyakov loop 
and conventional chiral condensate increases as we in- 
crease the current quark mass (see Fig. [7]). Initially the 
difference is zero for zero current quark mass but for 
m = 200MeV we find about 26MeV difference between 
the two temperatures. It is worthy of mentioning that 
this result is just for illustrative purpose as there are 
limitation of using NJL model with such a huge current 
quark mass. 
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and fi = at three temperature cases: 

T = T c x and T > T c x , where T c x is the chiral 



100 
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around cf> — n increases monotonously as temperature 
increases. As a result, dSi/dT increases. Above T c x , 
d{a)cf,/dT in the center region becomes zero (or very small 
in the case with finite current quark mass) and the region 
with large values of d{a) < f 1 /dT shrinks. Therefore, the 
integral in Eq. §5§ i.e. dEi/dT decreases as temperature 
increases. In all, dSi/dT gets its maximum at T* and 
the two transition temperatures T® and T* coincide in 
the case of second order chiral phase transition. 



FIG. 7: The critical temperatures T x and T c for different 
current masses. 



IV. THE RELATION BETWEEN T x AND Tf 

In the following, we offer a possible understanding on 
the simultaneity of the transition temperatures for 1st 
and 2nd order chiral phase transitions and the apparent 
difference between the two for the case of crossover. 

As mentioned earlier the transition temperatures are 
determined from the slope analysis of the conventional 
chiral condensate and the dressed Polyakov loop. So let 
us look at the temperature derivative of the general chiral 
condensate d(a}$/dT as functions of </>, the integral on 
which gives the temperature derivative of the Dressed 
Polyakov loop 



~dT 



# -id>d(o-)4, 



2tt 



dT 



(5) 



Fig. [8] shows d((j)tp/dT at different temperatures in the 
cases of second order chiral phase transitions. Around 
T^, large values of d{a) c f,/dT appear around <f> = it and 
dominate the integral in Eq. §5§. Below Tj? , d{o)<j,/dT 
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FIG. 9: Temperature derivative of the general chiral conden- 



sate 



dT 



for m = and /j, = 300MeV at three temperature 
cases: T <T C X ,T = T x and T >T C X , where T x is the chiral 
transition temperature. 



In the case of first order chiral phase transition, the 
situation is more complicated. Due to the discontinuity 
of {<7}<p at the first order chiral phase transition point, 



the temperature derivative of the dressed Polyakov loop 
can be expressed as 



~dT 



^-i<j> d ( a )>f> 



Cose 



2tt 



dT 



/c dT 1 



(6) 



where, the first term is determined by the regular behav- 
ior of d(a} f p/dT (see Fig^J), the second term is due to 
A(er) c , the jump of {a)^ at the first order phase transi- 
tion point at <j) = <fi c . When T < T x , the second term 
vanishes. Now we consider two limiting cases. First, in 
the case of a weak first order phase transition, A(cr) c is 
small and d{o) $ / dT around <j> = 7r is large, as showed in 
FiglU So the first term dominates the result of Eq. © 
and gives the similar result as that in the case of a sec- 
ond order chiral phase transition. Second, in the case of 
a strong first order chiral phase transition, A(c) c is large 
and d{a)tj>/dT is small. So the second term dominates 
the result of Eq. (JSJ. Then =jh. jg strongly dependent 
on the detailed information of d<j) c /dT. Our numerical 
results show that d(f> c /dT decreases as temperature in- 
creases, and so the second term also gives a decreasing 
contribution. In all, it is clear that dY>i/dT gets its max- 
imum at T x , i.e. T<? and T x coincide in the case of a 
weak first order chiral phase transition due to remnants 
of second order chiral phase transition. The coincidence 
in the case of a general first order chiral phase transi- 
tion is supported by our numerical results and can be 
generally expected. 
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FIG. 10: Second derivative of the general chiral condensate 

J^.^ for m = 5.5MeV and /x = at three temperature 
cases: T < T c x , T = T c x and T > T c x , where T c x is the 
chiral transition temperature. 



For the case with finite quark mass and small chemical 
potential, we have no phase transitions but crossover. So 
let us consider the second temperature derivative of the 
dressed Polyakov loop 



d 2 Si 
dT 2 



-i<p 



2tt 



d 2 (a)$ 
dT 2 : 



(7) 



whose zero point determines the transition temperature 
Tjf. Fig. [TU] shows the second derivatives of the gen- 
eral chiral condensate d 2 (a) <f, / dT 2 at three temperature 
cases: T <T C X ,T = T c x and T >T C X , where T c x is the 
chiral transition temperature. Similar to our previous 
observation, large values of d 2 (a) $ j dT 2 appear around 
4> = ir (see the two maximums in Fig. 1 10[) . as remnants 
of the second order phase transition in chiral limit. The 
difference is that d 2 (a) l j > /dT 2 in the center region (see 
the minimum in Fig. I10[) , is suppressed below T c x , ap- 
proaches zero at T c x and changes its sign above T c x . For 
T < T c x , d 2 (a) r p/dT 2 around the two maximums domi- 
nate the integral in Eq.© and d 2 T,i/dT 2 does not change 
its sign. Above T c x , the negative part around the mini- 
mum cancels the contributions from the maximums, and 
up to a certain temperature Tj 3 , this cancelation leads to 
the zero of the integral in Eq.Q. In all, the zero point 
of d 2 Yji/dT 2 comes from the negative contribution of the 
minimum at T > T c x , so the transition temperature re- 
lated with the dressed Polyakov loop must be greater 
than the chiral transition temperature, i.e. T c > T c x . 



V. CONCLUSION AND DISCUSSION 

We investigate the chiral condensate and the dressed 
Polyakov loop or dual chiral condensate at finite temper- 
ature and density in the two-flavor Nambu-Jona-Lasinio 
model. We find the behavior of dressed Polyakov loop in 
absence of any confinement mechanism still shows an or- 
der parameter like behavior. It is found that in the chiral 
limit, the critical temperature for chiral phase transition 
coincides with that of the dressed Polyakov loop. In the 
case of explicit chiral symmetry breaking, it is found that 
the critical temperature for chiral transition T x is smaller 
than that of the dressed Polyakov loop T® in the low 
baryon density region where the transition is a crossover. 
With the increase of current quark mass the difference 
between the two critical temperatures is found to be in- 
creasing. However, the two critical temperatures coincide 
in the high baryon density region where the phase tran- 
sition is of first order. 

From symmetry analysis, the dressed Polyakov loop 
can be regarded as an equivalent order parameter of de- 
confinement phase transition for confining theory. In 
the NJL model, the gluon dynamics is encoded in a 
static coupling constant for four point contact interac- 
tion. Since in this work we have included only quark 
degrees of freedom, a quantitative comparison will not 
match with other results. But the interesting fact is 
that the qualitative features (angle variation, tempera- 
ture variation) of the dressed Polyakov loop remains the 
same. Moreover, we expect that independent of the in- 
put of gluedynamics to the quark propagator, it is a gen- 
eral feature for T x = Tjf in the case of 1st and 2nd 
order phase transitions, and T x < T^P in the case of 
crossover, which qualitatively agrees with the lattice re- 
sult in Ref. [H|. This might indicate that for full QCD, 



in the crossover case, there exists a small region where 
chiral symmetry is restored but the color degrees of free- 
dom are still confined. This result should be checked in 
other effective models, e.g. in the framework of Dyson- 
Schwinger equations (DSE). 

The (T, fx) phase diagram with three flavors and with 
Ua(1) anomaly will be studied in the near future. 
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